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We consider Higgs boson production through gluon-gluon fusion in hadron colhsions. We 
present a numerical program that computes the cross section up to NNLO in QCD pertur- 
bation theory. The program aUows the user to apply arbitrary cuts on the momenta of the 
partons and of the photons or leptons that are produced in the final state. We present selected 
numerical results at the Tevatron and the LHC. 



1 Introduction 

The Higgs boson is a fundamental ingredient in the Standard Model (SM) but it has not been 
observed yet. Its search is currently being carried out at the Tevatron, and will be continued 
at the LHC. Once the Higgs boson is found, the LHC will be able to study its properties like 
couplings and decay widths. 

At hadron colliders, the SM Higgs boson is mainly produced by gluon-gluon fusion through 
a heavy-quark loop. For such an important process, it is essential to have reliable theoretical 
predictions for the cross section and the associated distributions. At leading order (LO) in 
QCD perturbation theory, the cross section is proportional to ag, as being the QCD coupling. 
The QCD radiative corrections to the total cross section are known at the next-to-leading order 
(NLO)^ and at the next-to-next-to-leading order (NNLO)^. The effects of a jet veto on the 
total cross section has been studied up to NNLO^^. We recall that all the results at NNLO have 
been obtained by using the large- approximation, Mt being the mass of the top quark. 

These NNLO calculations can be supplemented with soft-gluon resummed calculations at 
next-to-next-to-leading logarithmic (NNLL) accuracy either to improve the quantitative accu- 
racy of the perturbative predictions (as in the case of the total cross sectionSEl) or to provide 
reliable predictions in phase-space regions where fixed-order calculations are known to failEI. 

These fixed order and resummed calculations share a common feature: they are fully inclusive 
over the produced final state (in particular, over final-state QCD radiation). Therefore they refer 



to situations where the experimental cuts are either ignored (as in the case of the total cross 
section) or taken into account only in simplified cases (as in the case of the jet vetoed cross 
section). The impact of higher-order corrections may be strongly dependent on the details of 
the applied cuts and also the shape of various distributions is typically affected by these details. 

The first NNLO calculation that fully takes into account experimental cuts was reported in 
Ref.^, considering the decay mode H 77. In Ref.'Sj the calculation is extended to the decay 
mode H WW lulu (see also^^). 

In Ref.^I^ we have proposed a method to perform NNLO calculations and we have applied 
it to perform an independent computation of the Higgs production cross section. The calcu- 
lation is implemented in a fully-exclusive parton level event generator. This feature makes it 
particularly suitable for practical applications to the computation of distributions in the form 
of bin histograms. Our numerical program can be downloaded from The decay modes that 
are currently implemented are H ^ ^^^^ , — > WW — > lulu and H ZZ — > 4 leptons^I2]_ 

In the following we present a brief selection of results that can be obtained by our program. 
We use the MRST2004 parton distributions^^, with parton densities and as evaluated at each 
corresponding order (i.e., we use (n -|- l)-loop as at N'^LO, with n = 0,1,2). Unless stated 
otherwise, the renormalization and factorization scales are fixed to the value = ^xp = Mh, 
where Mh is the mass of the Higgs boson. 

2 Results at the LHC 

We consider the production of a Higgs boson with mass Mh = 165 GeV at the LHC {pp collisions 
at ^ = 14 TeV) in the decay mode H WW lulu The NLO and NNLO inclusive K- 
factors are Knlq = 1-84 and Kp^]\fLo = 2.21, respectively. We apply a set of selection cuts 
taken from the study of Ref.^^. The charged leptons (with \rj\ < 2) are classified according to 
their maximum (pti) and minimum {pt2) transverse momentum. Then pT2 should be larger 
than 25 GeV, and pTi should be between 35 and 50 GeV. The missing pT of the event is required 
to be larger than 20 GeV and the invariant mass of the charged leptons must be smaller than 35 
GeV. The azimuthal separation of the charged leptons in the transverse plane (At/)) is smaller 
than 45°. Finally, there should be no jet with pj}" larger than 30 GeV. 
The corresponding cross sections are reported in Table [TJ 



a{fh) 


LO 


NLO 


NNLO 


tJ'F = tJ'R = Mh/2 


17.36 ± 0.02 


18.11 ±0.08 


15.70 ±0.32 


fJ-F = tJ'R = Mh 


14.39 ± 0.02 


17.07 ± 0.06 


15.99 ±0.23 


llF = liR = '^Mh 


12.00 ± 0.02 


15.94 ±0.05 


15.68 ±0.20 



Table 1: Cross sections forpp H + X ^ WW + X ^ Ivlv + X at the LHC when selection cuts are applied. 

The cuts are quite hard, the efficiency being 8% at NLO and 6% at NNLO. The scale 
dependence of the result is strongly reduced at NNLO, being of the order of the error from the 
numerical integration. The impact of higher order corrections is also drastically changed. The 
X-factor is now 1.19 at NLO and 1.11 at NNLO. As expected, the jet veto tends to stabilize 
the perturbative expansion, and the NNLO cross section turns out to be smaller than the NLO 
one. The study of Ref. ^ shows that the efficiencies obtained at NNLO are in good agreement 
with those predicted by the MC@NLO event generator^^^. 

We now consider the production of a Higgs boson with mass Mh = 200 GeV'-^^'. In this 
mass region the dominant decay mode is — > ZZ — > 4 leptons, providing a clean four lepton 
signature. In the following we consider the decay of the Higgs boson in two identical lepton 
pairs. When no cuts are applied, the NLO i^T-factor is Knlq = 1-87 whereas at NNLO we have 



Knnlo = 2.26. 

We consider the following cuts'^. The transverse momenta of the leptons, ordered from 
the largest {pxi) to the smallest (pta) are required to fulfil pTi > 30 GeV, pT2 > 25 GeV, 
PT3 > 15 GeV and pxi > 7 GeV. The leptons should be central {\ri\ < 2.5) and isolated: the 
total transverse energy Et in a cone of radius 0.2 around each lepton should fulfil Et < 0.05 pT- 
For each possible e^e~ pair, the closest (mi) and next-to-closest {1112) to Mz are found. Then 
mi and m2 are required to be 81 GeV < mi < 101 GeV and 40 GeV < m2 < 110 GeV. The 
corresponding cross sections are reported in Table [2l 



a{fh) 


LO 


NLO 


NNLO 


^■F = l^R = Mh/2 


1.541 ± 0.002 


2.764 ± 0.005 


3.013 ±0.023 


fJ-F = tJ'R = Mh 


1.264 ±0.001 


2.360 ± 0.003 


2.805 ±0.015 




1.047 ±0.001 


2.044 ± 0.003 


2.585 ±0.010 



Table 2: Cross sections for pp ^ H + X ZZ + X — ^ e e e+e + X at the LHC when cuts are applied. 

Contrary to what happens in the H WW — > li^li^ decay mode, the cuts are quite mild, 
the efficiency being 63% at NLO and 62% at NNLO. The NLO and NNLO ivT-factors are 1.87 
and 2.22, respectively. Comparing with the inclusive case, we conclude that these cuts do not 
change significantly the impact of QCD radiative corrections. 
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Figure 1: Tranverse momentum spectra of the final state leptons for pp ^ H + X ~* ZZ + X ^ e^e~e^e~ + X, 
ordered according to decreasing pr , at LO (dotted), NLO (dashed), NNLO (solid). 

In Fig. [1] we report the pT spectra of the charged leptons. We note that at LO, without cuts, 
the Pti and pT2 are kinematically bounded by Mh/2, whereas pT3 < Mh/3 and pxi < Mh/4:- 
As is well known, in the presence of a kinematical boundary, fixed order predictions may develop 
perturbative instabilities^!^ at higher orders. In the present case, the distributions smoothly 
reach the kinematical boundary and no perturbative instability is observed beyond LO. 

3 Results at the Tevatron 

We now consider the production of a Higgs boson of mass Mjj = 160 GeV at the Tevatron (pp 
collisions at ^/s = 1.96 TeV). We consider the decay mode H — > WW — > lul'i^' with 1,1' = e,/Li 



and use the cuts from^^^. The inclusive JC-factors at NLO and NNLO are K^lo = 2.42 and 
Knnlo = 3.31, respectively. The trigger requires either a central electron with |?7| < 1.1 and 
Et > 18 GeV, or a forward electron with 1.2 < |Ty| < 2 with Et > 20 GeV and > 15 GeV, 
or a central muon with |7/| < 1 and px > 18 GeV. The leptons should be isolated: the energy in 
a cone of radius R = 0.4 around each lepton should fulfill E < 0.1 px- The selection cuts require 
PTi > 30 GeV, PT2 > 25 GeV and > 40 GeV. If < 50 GeV the azimuthal separation 
A(f){^T,p) should be larger than 20° for each lepton or jet. The invariant mass of the charged 
leptons should be between 16 and 75 GeV. The scalar sum of the pT of the leptons and of 
should be smaller than Mh- Finally, jets with E^ > 15 GeV and \7]\ < 2.5 are counted. We 
require either no such jet, or one of such jets with Et smaller than 55 GeV, or two of such jets 
with Et smaller than 40 GeV. Since the small A(j) region is the most important, we further 
require Acp < 80°. 

With these cuts the LO, NLO and NNLO cross sections are aio = 1-571 ib 0.003, onlo = 
3.16 ± 0.01, (Jnnio = 2.78 ± 0.17 fb. The corresponding ivT-factors are K^lq = 2.01 and 
Knnlo = 1.77, and exibit a strong reduction with respect to the inclusive case. As for the LHC, 
the selection cuts drastically change the impact of radiative corrections. Also the efficiencies are 
reduced: comparing the accepted and inclusive cross sections we obtain: e^o = 33%, eNLO = 
27% and e^NLO = 18%. These results suggest the existence of large theoretical uncertainties 
that need to be further investigated. 
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